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The three-dimensional distribution of astronomical objects observed in redshift space 
significantly differs from the true distribution, since the distance to each object cannot be 
determined by its redshift z only; for z< 1 the peculiar velocity field contaminates the true 
recession velocity of the Hubble flow, while the true distance for objects with z ^ 1 sensitively 
depends on the (unknown and thus assumed) cosmological parameters. This hampers the 
effort to understand the true distribution of the large-scale structure of the universe. In 
ON ■ addition, all cosmological observations are carried out on a light-cone, the null hypersurfacc 

of an observer at 2 = 0. This implies that their intrinsic properties and clustering statistics 
should change even within the survey volume. Therefore, a proper comparison taking account 
of the light-cone effect is important to extract any cosmological information from redshift 
catalogues, especially for z £ 1. We present recent theoretical developments on the two 
effects - the cosmological light-cone effect and the cosmological redshift-space distortion - 
that should play key roles in observational cosmology in the 21st century. 



§1. Introduction 



Galaxy redshift surveys in 1980s revealed and established the existence of large- 
scale structured extending around ~ lOOMpc in the current universe at z = 0. 
Theoretically, many cosmological models are known to be more or less successful in 
reproducing the structure at redshift z ~ 0. In fact, however, this may be largely 
because there are still several degrees of freedom or cosmological parameters needed 
to appropriately fit the observations at z ~ 0, including the density parameter, Qq, 
the mass fluctuation amplitude at the top-hat window radius of 8/i _1 Mpc, as, the 
Hubble constant in units of lOOkm/sec/Mpc, h, and even the cosmological constant 
Ao- This kind of degeneracy in cosmological parameters among viable models can be 
broken by combining the data at higher z. 

With the on-going redshift surveys of millions of galaxies and quasars and with 
large telescopes with high spectral resolution, one can probe directly the epoch of 
galaxy formation. One of the most important goals of cosmology in the next century 
is to construct a physical model of galaxy formation and evolution in the observa- 
tionally determined cosmological context. To this time this process has been simply 
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parameterized by the notorious bias parameter b, whatever its meaning might be. 
Presently, many theoretical and observational attempts are in progress to replace 
the parameter b by another physical model. Naturally, observational explorations of 
the larger-scale structure at z = and higher redshifts provide important clues to 
understanding the origin of structure in the universe. 

Redshift surveys of galaxies definitely serve as the central database for obser- 
vational cosmology. In addition to the existing catalogues including CfAl, CfA2, 
SSRS, and the Las Campanas survey, upcoming surveys such as 2dF and SDSS are 
expected to provide important clues to our universe. In addition to those shallower 
surveys, clustering in the universe in the range z = 1 — 3 has been partially revealed 
by, for instance, the Lyman-break galaxiesEP and X-ray selected AGNs.lf In partic- 
ular, the 2dFi§ (2-degree Field Survey) and SDSS (Sloan Digital Sky Survey) QSO 
redshift surveys promise to extend the observable scale of the universe by an order 
of magnitude, up to a few Gpc. A proper interpretation of such redshift surveys in 
terms of the clustering evolution, however, requires an understanding of many cos- 
mological effects which can be neglected and thus have not been considered seriously 
in redshift surveys of z <C 1 objects. 

This paper consists of two topics which should play key roles in the theoretical 
interpretation of the future redshift surveys of high-redshift objects, the cosmological 
light-cone effect (§^j) and redshift-space distortion (§||). Primarily, we intend to 
review and describe the two effects in a systematic and comprehensive manner on 
the basis of several of our papers. EP ~ Ef In addition, however, we input new materials 
in §3.4 and §3.5. Also, §2.2 presents theoretical predictions based on a different bias 
model from that adopted in one of our previous studies. & 

In this spirit, the remainder of this paper is organized as follows. Section 2.1 
briefly outlines a theoretical formulation of the two-point correlation function on 
the light-cone hypersurface, following Ref. ||). The corresponding theoretical predic- 
tions are presented in Section 2.2 with future QSO redshift surveys in mind. The 
predictions are based on a different model for evolution of bias from that adopted 
in Ref. ||). Thus they illustrate the extent to which the effect of bias changes the 
observable clustering of high-redshift objects. Section 2.3 summarizes the light-cone 
effect on the higher-order clustering statistics following Ref. ^) . Section || starts with 
the basic idea of the cosmological redshift-space distortion (§3.1) and its formulation 
in linear theory, both of which are on the basis of Ref. 0) . Then we comment on the 
systematic bias in estimating the cosmological parameter from shallower (z 0.2) 
galaxy redshift surveys, following Ref. |8|). The next two subsections are entirely new; 
Section 3.4 considers uncertainties due to the distance formulae in inhomogeneous 
cosmological models and to the evolution model of bias, and Section 3.5 examines 
the feasibility of the cosmological redshift-space distortion as a cosmological test 
to probe Qq and Aq- In the latter we fully explore the nonlinear effects also using 
high-resolution iV-body simulations. Finally, we summarize the main conclusions in 
Section ||. 
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§2. Cosmological light-cone effect 



Observing a distant patch of the universe is equivalent to observing the past. 
Due to the finite light velocity, a line-of-sight direction of a redshift survey is along 
the time, as well as spatial, coordinate axis. Therefore the entire sample does not 
consist of objects on a constant-time hypersurface, but rather on a light-cone, i.e., a 
null hypersurface defined by observers at z = 0. This implies that many properties of 
the objects change across the depth of the survey volume, including the mean density, 
the amplitude of spatial clustering of dark matter, the bias of luminous objects with 
respect to mass, and the intrinsic evolution of the absolute magnitude and spectral 
energy distribution. These aspects should be properly taken into account in order 
to extract cosmological information from observed samples of redshift surveys. 

For the CfA galaxy survey, for instance, the survey depth extends up to a 
recession velocity of 15000 km/s, which is interpreted as either <i max = 150h Mpc 
in spatial distance or z max = 0.05 in time difference. This translates to a ~ 10% 
difference in the amplitude of £ and P(k) in linear theory. Compared with the 
statistical error of the available sample, this level of systematic effect is negligible. 
Thus it is quite common to compare the observed £ with the theoretical predictions 
at z = 0. The situation will be entirely different for the upcoming galaxy and QSO 
redshift surveys, 2dF and SDSS; 0.3 — 1 million galaxies up to z max = 0.2, and 
0.3 — 1 x 10 5 QSOs up to z max = 3 — 5. Such observational samples motivate us 
to formulate a theory to describe the clustering statistics, fully incorporating the 
light-cone effect. In the remainder of this section, we present theoretical predictions 
for two-point & and higher-order correlation & functions which are properly defined 
on the light-cone. 

2.1. Defining two-point correlation functions on a light-cone 

In this subsection, we derive an expression for the two-point correlation function 
on the light-cone hypersurface in the spatially-flat Friedmann - Robertson - Walker 
space-time for simplicity; the line element is given in terms of the conformal time 77 

as 



Since our fiducial observer is located at the origin of the coordinates (rj = tjq, r = 0), 
an object at r and rj on the light-cone hypersurface satisfies the simple relation 
r = 7] -r]. 

We denote the comoving number density of observed objects (galaxies or QSOs 
satisfying the selection criteria) at rj and x = (r, 7) by 71(77, x )- Then the correspond- 
ing number density defined on the light-cone is written as 



If we introduce the mean observed number density (comoving) and the density fluc- 
tuation at rj, 710(77) and A(r], x), on the constant-time hypersurface, 



ds 2 = a 2 {r/) -drj 2 + dr 2 + r 2 dQ 2 . 



(2-1) 



77 LO (r,7) 



77,(770 - r, r, 



7)- 



(2-2) 



71(77, x ) = n o(v) I 1 + A(r), x)] , 



(2-3) 
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Eq. ( |2-2| ) can be rewritten as 

n LC (r, 7) = ra (% - r) [1 + ^(vo ~ r, r, 7)] . (2-4) 

The observed number density 710(77) is different from the true density of the objects 
71(77) at 77 by a factor of the selection function <f>(rj): 

n (r/) =n(ri)(/)(rf). (2-5) 

Thus no (77) already includes the selection criteria, which depend on the luminosity 
function of the objects and thus the magnitude-limit of the survey, for instance: 
When n LC (r, 7) is given, one may compute the following two-point statistics: 

X{R) = ^Lc / ^ / rldndQ^ J rldr 2 dQ, n 

x n LC (r 1 ,7 1 )n LC (r 2 ,7 2 )^(x 1 - x 2 - R). (2-6) 

Here xi = (ri, 7*17*1), x 2 = (r 2 , r 2 7 2 ), R = |R|, R = H/R, and V LC is the comoving 
survey volume of the data catalogue: 

V = / r 2 dr / dfy = — (r max 3 - r min 3 ), (2-7) 

with r max = r(z maiX ) and r min = r(z m - m ) the boundaries of the survey volume. Al- 
though the second equality as well as the analysis below assumes that the survey 
volume extends to 4tt steradian, all the results below can be easily generalized to the 
case of the finite angular extent. 

Substituting Eq. (2-4), the ensemble average of an estimator X{R) can be ex- 
plicitly written as 

(X{R)) = W(R)+U(R), (2-8) 

where 

W{R) = y^c / J drir i / d ^~n J dr ^ r l J d% 2 n (y7o - ri)n (?7o - r 2 ) 
x /z\(77o - ri, n, 7i)Z\(r/ - r 2 , r 2> 72)^^ (x x - x 2 - R), (2-9) 



and 



x 710(770 - rx)n (rjo - r 2 )<5 (3) (xi - x 2 - R). (2-10) 

After a tedious but straightforward calculation, B we have shown that the above 
definitions can be approximated as 

47T fmax 

W ( jR ) - TTlc / r dr M% - r)] 2 £(R; 770 - r) S ourcc , (2-lla) 

* » Tin in 



m 1 n 
max 
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where (,(R; ?/)source is the conventional two-point correlation defined on the constant 
hypersurface at the source's position. Note that U (R) is independent of R for R -C 
r max , as expected. 

The next task is to define the two-point correlation function on the light-cone. 
We propose the definition 



, P , , drr n (i] - r) £(R; r/ - r) SoU rce 

Z LC (R) = j2§ = ^ 7*= , (2-12) 

U ^ R > / drr 2 n (?7 - rf 

J ^min 

where the second expression uses Eqs. (2-lla) and (2-llb). If the correlation function 
of objects does not evolve, i.e., £(R; r] - r)s urce = (,(R; ??o)source, Eq. ( |2 • 1 2|) readily 
yields 

£ LC (fl)=£OR;%)Source, (2-13) 

as should be the case. 

Equation (^T3j) can be directly evaluated from any observed sample. First, 
average over the angular distribution and estimate the differential redshift number 
count dN/dz of the objects. Second, distribute random particles over the whole 
sample volume so that they obey the same dN/dz. Then the conventional pair- 
count between the objects and random particles yields X{R) (although not (^X(R)^, 
of course), while U(R) can be estimated from the pair-count of the random particles 
themselves. 

2.2. Predicting two-point correlation functions on a light- cone 

The corresponding theoretical predictions can be easily computed also, once a 
set of cosmological parameters and a model for the evolution of bias are specified. 
To illustrate the behavior of the two-point correlation functions on the light-cone, 
we adopt the following models. 

(i) cosmological parameters: We consider three models based on cold dark mat- 
ter (CDM) cosmogonies; SCDM with (J? , A Q) h, a s ) = (1.0,0.0,0.5,0.56), and 
LCDM with (f2o, Ao, h, o~%) = (0.3,0.7,0.7,1.0). The normalization erg is deter- 
mined from the cluster abundances. & Then we use the following fitting formula© 
for the linear power-spectrum of mass fluctuation: 

AIM - AHl + 2.U q )/2Ma? 

lV ; y/1 + 3.89g + plgp + (5.46g) s + (6.71g) 4 V ' 

q = ^/(r/iMpc- 1 ), (2-14b) 

r = J? /i(T 7 o/2.7 K)- 2 exp[-^ b (l + V2h^ 1 )]. (2-14c) 

(ii) mass correlation function: Gravitational nonlinear evolution is included by 
using the following fitting formulae IHP'EJ) for the mass power spectrum: 

u 

(2-15a) 



[l + Al L (k NL )} 



1/3 ' 
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Jnl{x) 

A 
B 

a 

V 
V 

n eS (k) 



f~N L [A 2 L (k L )}, 



1 + B(3x + (Ax 



n[(ii)V(z)/W 

0.482(1 + n/3)" - 947 
0.226(1 + ra/3)~ L778 
3.310(1 + n/3)-° 244 
0.862(1 + n/3)~ - 287 
11.55(1 + n/3)~ - 423 
d\uA 2 L {k) 



din k 



3. 



k=k L 



Then the nonlinear mass correlation function is computed via 







00 .9 . sinkx dk 



(245b) 

(245c) 

(245d) 
(245e) 
(245f) 
(245g) 
(245h) 

(245i) 



(246) 



(iii) evolution of bias: This is by far the most_uncertain factor in the current 
modeling. We simply use a linear bias modelli 3 ! on the basis of perturbation 
theory: 

6(^ = 1 + ^1(60-1), (247) 



D(z) 



where bo is the present value of the bias parameter. We denote by D(z) the linear 
growth rate (normalized as D(z) = 1/(1 + z) for z — > oo): 



D(z) 



1 + z' 



H(z') 3 



dz'. 



Here H(z) is the Hubble parameter at redshift z: 



(248) 



H{z) = F V^o(l + zf + (1 - fib - A )(l + zY + A . (249) 

According to this simplified scheme, z) Source is given by b(z)S >NL (x; z). 

(iv) selection function: With a magnitudedimitedQSO sample in mind, we adopt 
the following Bd^and quasar luminosity function. IliP '23 For 0.3 < z < 3, 

<P(M B ,z) = 1o0 .4 {q+1)[Mb _ m , (z)] - 10 o.4(/3+l)[M B -M*( 2 )] ' O 20a ) 

M£(z) = M^-2.5K L log(l + z), (2-20b) 



with = -20.91 + 51og/i, k l = 3.15, a 
10- 6 /i 3 Mpc 3 . For z > 3, 



-3.79, /? 



■1.44, <2>* 



6.4 x 



^( M B,z) 10 0.4(a+l)[M B -M*( 2 )] + 10 0.4(^+1)[A/ B -M*( Z )] ' ( 2 ' 21a ) 

Mg(z) = - 2.5K L log4 + S, (2-21b) 
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with A = (z — 3) log 3.2, B = 2.5A/(a — (3). To compute the B-band apparent 
magnitude from a quasar of absolute magnitude Mb at z (with the luminosity 
distance d L ), we apply the K-correction, 

B = M B + 51og(d L (z)/10pc) - 2.5(1 - p) log(l + z), (2-22) 

for the quasar energy spectrum L v oc v~ v with p = 0.5. While this luminosity 
function is derived from observed data assuming fi$ and Ao, we use this also for 
other cosmological models to make clean the differences due purely to the light- 
cone effect. 




10-' 1 10' 10 2 1 10' 10 a 

R [h-'Mpc] R [tr'Mpc] 



Fig. 1. Two-point correlation functions in the cluster normalized standard CDM model, (a) £ LC (R) 
denned on the light-cone hypersurface for QSOs with Bu m = 20. We assume three cases for 
the biasisc parameter, bo — 3, 2 and 1, from top to bottom. A nonlinear mass correlation 
functionli2Pis used for the solid lines, while the linear theory is used for dashed lines, (b) 
Linear (lower curves) and nonlinear Ex (upper curves) mass correlation functions defined on the 
constant-time hypersurfaces z = 0, 1 and 2. 




10-' 1 10' 10= 1 10' 10 z 

R [h-'Mpc] R [h-'Mpc] 



Fig. 2. Same as Fig. 1, but in a cluster normalized, low-density spatially flat CDM model. 
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The results are plotted in Figs. 1 and ^ for SCDM and LCDM models, re- 
spectively. The B-band limiting magnitude Bu m = 20 roughly corresponds to the 
upcoming SDSS QSO sample. The corresponding evolution of the amplitude of the 
correlation at R = 15/i _1 Mpc is plotted in Fig. 3. In the specific bias model we 
adopted, the amplitude monotonically decreases with increasing z. This is inconsis- 
tent with an observational claim that the QSO correlation amplitude increases 
as z. Given the theoretical uncertainties of the current theoretical understanding 
of the bias, it is premature to draw any decisive conclusion at this point. In fact, 
behavior consistent with the observational claim is obtained with a different model 
of bias. Nevertheless, this example illustrates the potential importance of the 
light-cone effect in understanding the evolution of the bias of high-redshift objects. 




0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 



z z 

Fig. 3. Evolution of amplitudes of two-point correlation functions at R = 15/i _1 Mpc of QSOs on 
the light-cone in cluster normalized (a) standard and (b) low-density spatially flat CDM models. 



2.3. Higher-order statistics on light-cone 

Let us move to the higher-order statistics of clustering on the light-cone hypersur- 
face. In particular we focus on the volume-averaged iV-point correlation functions, 
£n(R;z), at a redshift z and on a comoving smoothing scale R. In the higher- 
order statistics, it is more useful to introduce the normalized higher-order moments 
Sn(R; z) = £jv(i?; z )/[^2(R'i z )] 1 than £n(R; z). The hierarchical clustering ansatz 
states that Sn(R', z) is constant and independent of the scale R. Moreover, £n(R) z) 

evolves in proportion to|£ 2 (-R;2) > an d Sn(R; z) is independent of z, according 

to perturbation theory. S*'© 

As described in §2.1 and §2.2, however, the iV-point correlation functions aver- 
aged over the light-cone, 



z dz w(z)£ iN (R; z) 

~£n(R'i < z max) = — Jz^k ' (2-23) 

z 2 dz w(z) 
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and the corresponding moments, 
S N (R; < z B 



Zn(R; < 



^max, 



N-l 



(2-24) 



are the statistics more directly estimated from redshift surveys than their counter- 
parts defined on the idealistic z = hypersurface. Note that the expression ( 2-23| ) 



looks slightly different from Eq. ( 2-12 ). This is because we have count-in-cell analy- 
sis in mind with the sampling cells being placed randomly in the z-coordinate. The 
effect of the selection function is taken into account by the weight function w(z). In 
the case of count-in-cell analysis, one can correct for the selection function <f>(z) by 
multiplying the count in cells located at z by l/(p(z). Then w(z) can be set to unity 
for z < z max and zero for z > z max in principle, where z mSuX is the maximum redshift 
of the sample. In the remainder of this subsection, we assume that the effect of the 
selection function is already corrected in this way, and consider the light-cone effect 
on the moments ( |2-24 ) due to the difference of the gravitational evolution within the 
survey volume. 

Let us define a function G which describes the evolution of the averaged two- 
point correlation function at R and z: 



C(R;z) = G(R;z)^(R;0). 



(2-25) 



In general, G is not a simple function of R and z, but a complicated functional of 
£. In the linear regime, however, G is independent of R and given by [D(z)/ D(0)] 2 , 
and even in the nonlinear regime, it is known that G is approximately expressed 
as a function of R and z alone. To proceed more specifically, we apply the fitting 
formulae^ which relates the evolved two-point correlation function ^(-R; z) with its 
linear counterpart ^ l (Rq;z) as follows: 



F[x) 



B(n)F[Z L (Ro;z)/B(n)}, 
x + 0.45a; 2 - 0.02x 5 + 0.05z 6 
1 + 0.02x 3 + 0.003x 9 / 2 



(2-26a) 
(2-26b) 



In the above equations, n denotes the effective spectral index of the power spectrum 
evaluated at the scale just entering the nonlinear regime, Rq = [1 + £%(z, R^^R, 
and B{n) = [(3 + n)/3] - 8 . 

The inverse of Eq. (2-26a) is also empirically fitted as follows: 



f-\v) 



B(n)F- 1 [Z E (R;z)/B(n)], 

I 1 + 0.036y L93 + 0.0001?/ 3 



1/3 



y 



1 + 1.75t/ - 0.0015y 3 b3 + 0.028y 



(2-27a) 
(2-27b) 



Then the scale-dependent evolution factor G(z) 
expressed explicitly in terms of £#(i2; 0): 



G(R; z) defined by Eq. ( ^25|) is 



G(R;z)A E ^ Z 



B(n) 

l E {R;Q>) £ E (R;0y 



■F 



D 2 {z) p-i ( Ze(R,0) 
D 2 (0) { B(n) 



(2-28) 
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Substituting the evolution law (2-25), Eq. ( |2-24| ) is explicitly written as 



•3 m ax 



1 N-2 



Z max 



S N (R; < z max ) = I z z dz w(z) 

>o J ' Uo 

f ^max 

< / z 2 dz w(z)S N (R;z){G(R;z)} 
lo 



z 2 dz w{z)G(R;z) 



N-l 



N-l 



(2-29) 



In order to proceed further, we assume that Sn(R;z) does not evolve with z, 
i.e., Sn(R;z) = Sn(R;0). As described above, this is a reasonable approximation 
as long as objects are unbiased tracers of the underlying density field. Also let us 
introduce the measure of the light-cone effect: 



A N (R; < z n 



Sn(R~, < Z mi 

S N (R;0) 



l. 



(2-30) 



Then (1 + An) can be regarded as a correction factor for the Sn estimated without 
considering the light-cone effect. This quantifies the importance of the light cone ef- 
fect on the higher-order clustering statistics in future surveys. Using Eqs. (2-29) and 
( 2-28| ) and assuming Sn(z) = SV(0), we evaluate Aw(< for SCDM, LCDM, 

and OCDM, which have (fl , A , h, <x 8 ) = (1.0,0.0,0.5,0.56), (0.3,0.7,0.7,1.0), and 
(0.3, 0.7, 0.0, 1.0), respectively. Figure 4 displays Ajy(R; < z) as a function of z, while 
Figure 5 plots An(R; < z) against R. 



100 
10 

i 

o.i 

0.01 

100 
10 

1 

o.i 

0.01 

100 
10 

o. I 

0.01 



Y< <lllll — I uuj 




r OCDM ///A 
prVmni — 1 1 




- LCDM ///A 





0.1 
rcdshift z 



Fig. 4. log 10 An(R; z) are shown as functions 
of log 10 z at R = l/i _1 Mpc (left panels) and 
10/i _1 Mpc (right panels). The curves cor- 
respond to N = 3, . . . , N = 10 from bot- 
tom to top with iV = 3 and N — 7 plotted 
in thick lines (from Ref. H)). 




10 

Mpc) 



Fig. 5. log 10 An(R\ z) are displayed as func- 
tions of log 10 R at z = 0.2 (left panels) and 
1.0 (right panels). The curves correspond 
to N = 3, . . . , N = 10 from bottom to top 
with iV = 3 and N — 7 plotted in thick 
lines (from Ref. ^)). 



The figures suggest that the light-cone effect is quite robust. Although its details 
depend on the model, the difference is fairly small, and qualitatively all models 
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behave similarly; the magnitude of the correction monotonically increases for higher 
order N. Also as expected, the light-cone effect becomes larger as z max increases 
(Fig. 4). Although the correction is relatively small for shallow surveys with z ^, 0.2 
samples, An(R;< z max ) becomes ^ 10% in nonlinear scales (R ~ 1/t Mpc). In 
SCDM, for instance, A^(R; < z max ) exceeds unity for N > 6 for the entire dynamic 
range plotted. Furthermore, Fig. 5 indicates that even if the hierarchical ansatz 
is correct (i.e., Sn(R',z) is independent of R), the light-cone effect should generate 
apparent scale-dependence, since the correction behaves differently at different scales 
for a given redshift. In future surveys extending to z > 1, Fig. 5 implies that the 
required correction for the light-cone effect is appreciable, ranging from up to unity 
for S3 through factors of few for Sq to factors of hundred for Sio . 

§3. Cosmological redshift-space distortion 

3.1. Basic idea of cosmological redshift-space distortion 

The approach described in §|2] is based on an implicit idea to treat all objects 
in a survey catalogue simultaneously. If the number of objects in the catalogue is 
sufficiently large, one can divide the objects in many redshift bins. Then the light- 
cone effect discussed above is less important, as long as one treats each individual bin 
separately. In this case, however, another interesting effect due to the geometry of the 
universe emerges. This originates from the fact that the (observable) redshift-space 
separation s is mapped to the (unobservable) comoving separation x of objects at z 
differently, depending on whether the separation is parallel or perpendicular to the 
line-of-sight direction of an observer at z = 0. Due to this effect, a sphere located at 
z becomes elongated along the line-of-sight in general. c3 In this section, we describe 
the anisotropy in the two-point correlation function of high-redshift objects induced 
by this effect, which we call the cosmological redshift-space distortion, u>'c3 in order 
to dislingmsh the conventional redshift-space distortion due to the peculiar velocity 
field.!) -E) 

Throughout this section, we assume a standard Robertson - Walker metric of 
the form 

ds 2 = -dt 2 + a{t) 2 {d X 2 + S( X ) 2 [d9 2 + sin 2 Odcf 2 )}. (3-1) 

We adopt a normalization for which the present scale factor ao is unity. Then the 
spatial curvature K is related to the other parameters as 

K = H 2 (f2 + Ao - 1), (3-2) 

and S(x) is determined by the sign of K according to 

( sin (VKx)/VK (K>0) 
S(x) = { X (K = 0) . (3-3) 

{ smh(V^Kx)/V^K {K<0) 

The radial distance x( z ) is given by 
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Consider a pair of objects located at redshifts z\ and z-i- If both the redshift 
difference Sz = z\ — Z2 and the angular separation of the pair 59 are much less 
than unity, the comoving separations of the pair parallel and perpendicular to the 
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Fig. 6. Behavior of c\\{z) and c ± (z) for (f2 , A„) = (0.3,0.7), (1.0,0.0), and (0.3,0.0) corresponding 
to the solid, dashed and dotted lines, respectively. 
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Fig. 7. Effect of inhomogeneity on the angular diameter distance D A (z) (upper panels) and the 
correction factor n{z) = C\\ (z)/c±(z) (lower panels) for Ao = and 0.9 models in J?o = 1 (left 
panels) and J?o =0.1 (right panels) universes. Thick lines indicate the results for the empty 
beam (a = 0), while thin lines for the filled beam (a = 1). 
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line-of-sight direction, Xy and x ± , are given by 

*„(*) = ^* = S£, *x(*)=S(x(*))** = ^, (3-5) 

where c^(z) = Hq/H(z), c ± (z) = HqS(x( z ))/ z an d z = ( z l +^2)/2. Thus their ratio 
becomes 

X\\(z) C||(z) Sz _ . . Sz 

^7(zj = ^(z)zJe =1l{z) zJe- (3 ' 6) 

Since 5z/(z59) is the ratio of the parallel and perpendicular separations to the line- 
of-sight direction, rj(z) represents the distortion in the redshift space coordinates 
induced by the geometry of the universe. Typical behavior of c^(z), c ± (z) and rj{z) 
is plotted in Fig. 6 and the lower panels of Fig. 7. This is purely a general relativistic 
effect which was first proposed as a potential test of the cosmological constant. EIP 

Actually it took a couple of decades to find realistic phenomena to_wJhich the test 
could be applied observationally. Recently two independent groups lP'E3) proposed 
to use the anisotropy in the clustering pattern of quasars and galaxies at high- 
redshifts as a probe of r](z). The next subsections describe the methodology in 
linear theory, and then examine the feasibility in the nonlinear regime using iV-body 
simulations. IHZP 

3.2. Linear theory of cosmological redshift- space distortion 

We choose a fiducial point at redshift z as an origin, and set up locally Euclidean 
coordinates with respect to this point. Let us adopt the distant (plane-parallel) ob- 
server approximation and choose the line-of-sight direction as the third axis. Then, 
if an object is located at x = (xi,X2,xs) in the real (comoving) space, the corre- 
sponding redshift-space coordinates, s = (s±,S2,ss), observed at z = are written 
as 

si = X ) , (3-7a) 
s 2 = j-t, (3-7b) 



S3 



^obs * 



H c l{ (z) 



1 + z 
X3 + H{z) V " 



(3-7c) 



In the last expression, v\\ is the recession velocity of the object relative to the observer, 
and z <C H(z)x% is assumed. Computing the Jacobian of the above transformation in 
linear theory, one can relate the density contrasts of the objects in real and redshift 
spaces as 

«W( a (x)) = 5W(x)-^ftt; l . (3-8) 

The peculiar velocity in linear theory is written in terms of the mass density contrast 
<$ mass asH 

«,l(x) = -^£l/(z)a 3 A- 1 ( 5 mass (x), (3-9) 

1 ~~\~ Z 
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where A -1 is the inverse Laplacian, and 

f{z) = ,, ' ' - :)■■■" + 

n{z) = 



1 + 



Q{z) 



Hn 



i 2 



X(z) 



H(z) 

Hp 
H(z) 



(l + z) 3 f2 , 



1 2 



A 



o- 



(3- 10a) 
(340b) 
(3- 10c) 



In order to close the equations, one has to relate the density contrast of objects 
in real space 5^ to that of mass, <5 mass , by specifying the model of bias. As in §2.2, 
we adopt a linear bias: 

jW(x;z) = b(z)S mmm (x]z). (3-11) 



(3-12) 



Substituting the above equations into Eq. (|3-8j), we obtain 

d 3 k 



««(-(x)) 



(2tt) 2 



1 + 



e*' x (5W(k) J 



where /3(z) = f(z)/b(z), ancL(5^ r ^ is the Fourier transform of 5^ r \ Repeating the 
method of Hamilton (1992), c3) we obtain an explicit formula for the redshift-space 
two-point correlation function which is valid at arbitrary z in linear theory: 



z)] 2 )Ux)Po(p) 



|/3(z) + ^)] 2 ) 6(x)ft0*) 
+ ^[/3(^)] 2 C4(^)P4(M), 



1 

27 2 7 

[6(z)] 5 



dkk 2 j2i{kx)P(k; z), 

, n 2 



D(*) 



D(0) 



P mass (A;;z = 0). 



(3-13a) 
(3- 13b) 

(343c) 



In the above, x = ^ c^ 2 s^ 2 + c ± 2 s ± 2 , \i = C\\S\\jx (sy = S3 and s 2 
the P n are the Legendre polynomials, i.e., Pq{h) = 1, ^(m) = (3// 2 



a? + s|), 
l)/2, and 

P 4 ( A t) = (35^ 4 - 30/i 2 + 3)/8. Figure 8 plots £ (s) (s x ,S||) for T = 0.2 and cr 8 = 1.0 
models to illustrate the degree of distortion. Figure 9 shows the difference between 
CDM models with r of Eq. (244c) and fixed r{= 0.2) models with the same i?o 
and Ao- 

3.3. Implication for galaxy redshift surveys 

The cosmological distortion effect becomes important also even for shallower 
galaxy redshift surveys.© One may formally expand ^ s \s ± ,s n ) in terms of the 

observables, s = s 2 + s 2 ^ and fj, s = s^/s, instead of the unobservable variables 
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(x(z), fi x (z)): 



1=0 1=0 



(3-14) 



Since we are interested in surveys for which z <C 1, we can further expand the above 
summation up to first order in z, and then obtain 



^K^s^z) 



1=0 



The explicit expressions for Qi ( s 'i z ) °f the form: 



(^(x;z) = ^'(x;0)+S 2l z 



is), 



(3-15) 



(3-16) 



for 1 = 0, 1, 2 and 3 are found in Ref. ||). 

One possible application of those perturbative formulae is to estimate a system- 
atic error for the /3-parameter, /3q, due to the neglect of the cosmological redshift- 
space distortion. Hamilton (1992) proposed to estimate (3$ from the moments of the 
observed two-point correlation functions of galaxies on the basis of the relation© 



1 + 2/?q/3 + /? 2 /5 
4/3 /3 + 4/? 77 



^ s) (x;0) fx 
r 



3 dx 
x 



(3-17) 



lo e?(r;0) 

If one takes account of the cosmological redshift-space distortion at z, £q S \x; 0) and 
^ s (x; 0) on the right-hand side of the above equation should be replaced by (x; z) 
and C2 (^j z), respectively. Then substituting the perturbation expansion ( 3T6|) into 
Eq. ( |3T7D , one can compute the systematic error for {3q defined through 

P(z)=Po + ez. (3-18) 

The result i> consists of the two terms corresponding to the evolution of the (3- 
parameter and the geometrical effect: 



d/3(z) 



dz 



+ 



1 + . 



l + 6/? /7 + 3/3^/35 




1101 
15 



P 3 o 



1/4 + 3/3 /7 + 17/^/70 + $j/21 + $/196 9£(r; 0) 



/ r K(x;0)-^(r;0)]x2dx 



din', 



, (3-19) 



where qo is the deceleration parameter (= — \q). For magnitude-limited sam- 
ples, the above expression should be averaged over the sample with an appropriate 
weight according to the selection function. If the magnitude limit of the survey is 
18.5 (in B-band) as in the SDSS galaxy redshift sample, the systematic error ranges 
between —20% and 10%, depending on the cosmological parameters. Although such 
systematic errors are smaller than the statistical errors in the current surveys, they 
will definitely dominate the expected statistical error for future surveys. 
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3.4. Effects of the cosmological distance and evolution of bias 

Here we discuss two potentially important effects, which were ignored in §3.1 to 
§3.3, on the cosmological redshift-space distortion: the effect of inhomogeneities in 
the light propagation and the evolution of bias. 

The angular diameter distance D A , which plays a key role in the geometrical test 
at high z, depends sensitively on the inhomogeneous matter distribution as well as 
Ao and the density parameter Qq. A reasonably realistic approximation to the light 
propagation in an inhomogeneous universe is given by the Dyer - Roeder distance. E2P 
It assumes that the fraction a of the total matter in the universe is distributed 
smoothly and the rest is in clumps. If an observed beam of light propagates far from 
any clump, then the angular diameter distance D A {z; a, Q$, Ao) satisfies 



d 2 D A 
dz 2 + 



2 1 dH(z) 

T+z + H{z) dz 



dD A 3aH2{2 (l + z) 

+ o -ET^2 L d * = °> ( 3 - 2 °) 



dz 2 H{zf 

with D A (z = 0) = and dD A /dz(z = 0) = 1/Hq. The preceding discussion on 
the cosmological redshift-space distortion adopted a standard distance, which corre- 
sponds to the extreme case of a = 1. As shown in Fig. 7, the effect of inhomogeneity 
represented by the parameter a in the above approximation, however, could be large 
for high z if a is significantly different from unity. Another uncertainty will come 
from the possible time-dependence of the bias parameter b(z). As we emphasized 
in the discussion of the light-cone effect, we do not yet have any reliable theoretical 
model for bias. Thus we adopt the linear bias model fl2-17|) so as to highlight the 
effect of the evolution of bias in the present context. 

Figure 7 may seem to indicate that inhomogeneity makes even a larger difference 
than that of Ao, especially for z>l. In reality, however, the situation is not so bad; 
since the expectation value of a is determined by the effective volume of the beam 
of the light bundles, it depends on the depth z and the angular separation 59 (of 
the quasar pair in the present example). For larger z and larger 59, a(z,59) should 
approach unity in any case, and the result based on the standard distance as in §3.1 
to §3.3 would be more appropriate and closer to the truth. Since we do not have any 
justifiable model for a(z,59), we will consider two extreme cases, a{z,59) = 1 (filled 
beam) and (empty beam). Our main purpose here is to highlight the importance 
of the effect, even though more realistically a{z, 59) is somewhere between the two 
extreme cases; it is shown that a(z, 59) = 1 is a good approximation for z> 1 and 
59 > 1.0) 

It is quite reassuring that even in these extreme cases the inhomogeneity effect 
is much weaker than that of Ao up to z ^ 2 in low density universes, as the right 
panels in Fig. 7 illustrateclearly. Since a relatively low value of Qq around 0.1 — 0.3 is 
favored observationally,H the optimal redshift to determine Ao in low Q$ universes 
is z = 1 - 2. 



Figure 10 illustrated the evolution of bias (Eq. ( 2-17 ); upper panels) and of the 
resulting (3{z) parameter (lower panels). This implies that as long as Fry's model of 
b(z) is adopted, one can distinguish the value of Ao independently of the evolution 
of bias only in low density (Qq <C 1) models and at intermediate redshifts (z £,2). 
Together with the indication from Fig. 7, z = 1 — 2 would be an optimal regime to 
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probe Ao at least in low-density universes. Figure 11 illustrates the extent to which 
this is feasible simply on the basis of the anisotropy parameter, 

_ g('>(«n=o,*;*) r , 91 x 

£\s) ~&H8,8 ± = 0;z)' [6 ' 

adopting the power spectrum of the CDM models; in fi Q = 1 models the value of 
a completely changes the z-dependence of the anisotropy parameter while Qq = 0.1 
models are fairly insensitive to it. In addition, (s) / (s) f° r z ~ 2 in £2$ = 0.1 
models is basically determined by the biasing parameter at z = and less affected 
by the evolution of b(z). Figure 12 shows the scale-dependence of the anisotropy 
parameter in Qq = 0.1 and h = 0.7 CDM models. This clearly indicates that one 
can distinguish the different Ao and bias models by analyzing the anisotropy of the 
correlation function at z = 1 almost independently of a. 

3.5. Testing the redshift- space distortion with N-body simulations 

As illustrated in Figs. 7 to 9, the two-point correlation functions become elon- 
gated along the line-of-sight due to the cosmological redshift-space distortion in linear 
theory. In reality, the finger of god due to the non-linear peculiar velocity affects the 
distortion pattern in the same direction. Therefore the proper modeling of the non- 
linear effects is essential to estimate the cosmological parameters from the observed 
distortion. Also, the available number of observed objects would limit the statistical 
significance of the analysis. In order to examine these realistic effects in applying 
the redshift-space distor^mas a cosmological test, we use a series of high-resolution 
N-body simulations. e3)'E§)'EZP The simulations assume the three representative cos- 
mological models summarized in Table |. Each model has three realizations with 
different random seeds in generating the initial condition, and employs Np = 256 3 
dark matter particles in the simulation volume of (300/i _1 Mpc) 3 (comoving). Fig- 
ure 13 displays the results of ^ s ^(s ± ,S||) at z = 2.2 for these models. The upper 
panels plot the predictions in linear theory, the middle panels are computed from 
randomly sampled N = 5 x 10 5 particles, and the laffier_panels from the N = 2 x 10 4 
most massive halos (groups of particles) identified, E^'E^P so as to take into account 
the effects of the finite sampling and the biasing to some degree. 



Table I. Simulation model parameters. 



Model 


Oo 


Ao 


r 


0"8 


N 


realizations 


SCDM 


1.0 


0.0 


0.5 


0.6 


256 a 


3 


OCDM 


0.3 


0.0 


0.25 


1.0 


256 3 


3 


LCDM 


0.3 


0.7 


0.21 


1.0 


256 3 


3 



Figure 14 plots the reduced x 2 contours from ^ s \s ± , S\\) of simulations. Since 
our theoretical predictions do not include nonlinear effects at this point, we exclude 
the regions with s^/s± > 2, which are likely to be seriously contaminated by nonlinear 
peculiar velocities. While Fig. 14 demonstrates that the current methodology works 
in principle, the expected S/N is fairly low. This is largely because we adjusted the 
sampling rate for the high-z QSOs. The situation would be improved, though, if 
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we apply the present methodology to a statistical sample of Lyman-break galaxies, 
for instance, whose number density is larger and their strong clustering is already 
observed, a 1 

In order to examine the nonlinear effects in the cosmological redshift-space dis- 
tortion, we consider the power spectrum, rather than the two-point correlation func- 
tions, in which the phenomenological correction for the non- linear finger-of-god effect 
has already been discussed in the literature. i3' c3 Specifically, we model the power 
spectrum before the cosmological redshift-space distortion as 

P {s) (k,fj,) =P(k)[l + f3n 2 ] 2 D[kfMT P ], (3-22) 

where [x is the direction cosine in k space, and the second factor on the right-hand- 
side comes from Eq. ( |3-12| ). The last factor is a phenomenological correction for 
non-linear velocity effect. We assume that the pair- wise velocity distribution in real 
space is exponential with a constant pairwise peculiar velocity along the line-of-sight. 
<7p. In this case the damping term in Fourier space, D{kfio~p], is given byliiP 

Combining the geometrical effect, the power spectrum of objects at z observed 
in redshift space is expressed as 

p(CR ) ^fc s ^ k s u ; z) = t~To rT-^ (^J- > ^11 > z ) 
c±{z) 2 c ]] (z) 

= , V2 ( Z W + f3(z)v 2 } 2 D[k^p(.z)}, (3-24) 

where the relation of the comoving wave numbers in real space, k, and in cosmological 
redshift space, k s , is expressed as 

fen 



fc B = — rr, k = Jkl + k 2 , /i = -r, (3-25) 



c x (z)' c \\{ z )' ± k 

and P(k; z) is a comoving real-space power spectrum at redshift z. 

Clearly, the final expression for the redshift-space power spectrum p( CRD ) (k s± , k s ^ ; . 
depends on a number of parameters: J?o> \)> b(z), P(k;z), and crp(z). While 
none of these parameters has been determined precisely yet, there exist some tight 
constraints on them which can greatly reduce the number of the independent un- 
known parameters; provided that one adopts the linear biasing, the shape of the 
linear density power spectrum has already been determined fairly well by the APM 
galaxy survey, for instance. ® The upcoming redshift surveys of nearby galaxies will 
improve this measurement significantly. Then, given as and b(z), P(k;z) is already 
accurately determined. Furthermore, with future large surveys which are pertinent 
to the analysis here, it should be fairly easy to determine b(z) for a given cosmology. 
The pair-wise velocity dispersion o~p(R,z) at large separation can be determined by 
the other parameters through the cosmic energy equation:^ 

2 

o- p ,mjb = (vf 2 {r oo)) = g<>i>, (3-26a) 
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3Q(z)H 2 (z)I 2 (z) 
2(1 + z) 2 



dk A% L (k,z) 



1 + z 
D 2 (z) 



i 

1 + 2 



[l + z') 2 



dz 



k 



k 2 



(3-26b) 
(3-26c) 



Table II. Velocity dispersion ap from theoretical predictions and the present simulations (in units 
of km/sec). 



Model 


0\p,mjb(2 = 0) 


0-p,sim(z = 0) 


(Tp,ujb(z = 2.2) 


(Tp,sim(z = 2.2) 


SCDM 


580 


591 ±2 


164 


166 ±1 


OCDM 


599 


603 ±6 


368 


364 ±3 


LCDM 


593 


606 ± 10 


380 


377 ±4 



As Table II indicates, the above analytical fit is in good agreement with our 
simulation results. Finally, a constraint on erg and Qq from cluster abundances 
at z = is fairly well-established. Thus combining these model predictions and 
observational constraints, we will be left with only two unknown parameters, Qq and 
An, which we desire to determine from the cosmological redshift distortion. 



Figures 15 and 16 display the contour plots for p( CRD \k s± , k s ^) at z = and 
z = 2.2, respectively. The upper, middle, and lower panels correspond to theoretical 
predictions in linear theory, nonlinear model predictions on the basis of equations 
(2-15) and (3-24), and simulation results, respectively. Note that P(k) in this section 
is related to A 2 (k) in §2.2 as0> 

A 2 (k) = -L-47rk 3 P(k). (3-27) 

The middle panels plot two nonlinear models which adopt different ap in Eq. ( |3-23| ). 
The solid curves correspond to the pair-wise velocity dispersions directly evaluated 
from thesimulation data, while the dotted curves correspond to an analytical fitting 
formulaic (Eq. (3-26)). The right-hand-side of the above equation depends on the 
scale R through a spherical top-hat window function, W(k;R), while Eq. (3-23) is 



derived on the assumption that o~p(z) is scale-independent. Note that we adopt the 
velocity dispersion in comoving coordinates. This implies that we have to multiply 
the proper velocity by the conversion factor Hq(1 + z)/H(z). We adopt the value at 
R = 40/i~ 1 Mpc, which is the median value of the fitting range of our analysis (see 
below). 

As in the case of the two-point correlation functions, the degree to which one can 
recover the power spectrum sensitively depends on the number of sampled particles. 
The lower panels in Figs. 15 and[i~6| use all the simulation particles. We repeated the 
same analysis with randomly sampled Np = 5 x 10 5 , 5 x 10 4 and 5 x 10 3 particles, 
and the results are displayed in Fig. 17. The SDSS QSO surveys, for instance, expect 
to have O(10 4 ) QSOs between z = 1.5 and 2.5. This fi gure implies that although the 
phenomenological nonlinear models reproduce the simulation results very well, the 
statistical noise due to the limited numbers of QSOs will dominate the cosmological 
signal as long as one attempts to directly compare the p( CRD \k s± , k s \\)- 
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One can increase the signal-to-noise ratio by expanding the power spectrum in 
multipole moments as 

pf RD \k s ,z) ee *±I f^ CRD \k s ^ ]Z ) Pl {^ (3-28) 

where the Pi are the Legendre polynomials. In order to illustrate the higher signal- 
to-noise ratio in this approach, we plot the monopole term, PQ CRD \k s ; z = 2.2), in 
Fig. 18 computed from 5 x 10 4 randomly sampled particles. The quoted error bars 
are estimated from the la dispersions of P$ of 24 random subsamples in total (eight 
randomly sampled particle sets for three different realizations). The five curves of 
different line types correspond to theoretical predictions which use different values 
for ap quoted in the plot (but fix the other parameters as the values adopted in the 
simulations). While the power spectrum itself is rather noisy (see the middle panels 
in Fig. 17), the estimated moment is very robust. Figure 18 suggests that the best-fit 
ap is systematically smaller than the values listed in Table II. Better quantitative 



agreement is obtained by replacing the dp in Eq. ( 3-23 ) with the pair-wise velocity 
dispersion divided by y/2. This is related to the validity of the modeling of nonlinear 
velocity correction (e.g., Eq. ( 3-23Q ), and will be discussed in detail elsewhere. c3 



Finally we compute the x 2 between the theoretical predictions (with nonlinear 
corrections) and the simulations by varying the parameters. The results are shown as 
contours in Qq - ap and Qq - as planes in Figs. 19 and 20, respectively. In these fig- 
ures, the upper and lower panels correspond to the analysis based on p( CRD ) (fc s± , fc s y ) 

using all simulation particles (N = 256 3 ) and PQ RD \k s ) using 5 x 10 4 randomly 
sampled particles (one realization from each model in Table ||). The x 2_ fit is carried 
out in the range of (2vr/60)/iMpc _1 < k x ,k h k a < (27r/20)/iMpc _1 . 

In Fig. 19, we fix the value of as according to the fitting formula based on the 
following cluster abundances: EP 

as = (0.54 ± 0.02) x { ^ . 28 _ . 91fio+0 . 68fi , ^° " L ~ U ° h (3-29) 

"0 v^o — 0). 



In Fig. 20, we fix the value of ap according to Eq. (14) in Ref. [34]). Incidentally the 
cluster abundance constraints (3-29) are fairly orthogonal to our constraints from 
the redshift distortion. 

The best-fit values for ap and as in the above plots are slightly smaller than their 
true values (marked as crosses). This is related to the nonlinear velocity correction, 
as described above, and wecan easily correct for this systematic effect by adopting 
a more appropriate model. c3 Therefore we conclude that it is feasible to break the 
degeneracy in the cosmological parameters by combining the cosmological redshift- 
space distortion in the future QSO samples with other cosmological tests, despite 
the fact that the present modeling of nonlinear effects is fairly empirical. 
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§4. Conclusion 

The present paper focuses on two important effects, the cosmoiogicai light-cone 
effect and redshift-space distortion, which have been largely ignored in previous dis- 
cussions of clustering statistics. We have demonstrated that they play an important 
role in the analysis of the upcoming redshift surveys, particularly of high-redshift 
objects, as both cosmoiogicai signals and noise depending on specific aspects of the 
phenomena that one is interested in. We summarize our main conclusions here. 

1. We derived an expression for the two-point correlation function properly de- 
fined on the light-cone hyp er surf ace. |) This expression is easily evaluated nu- 
merically when the underlying cosmoiogicai model is specified. With this, one 
can directly confront the resulting predictions with the observational data in a 
fairly straightforward manner. 

2. The cosmoiogicai light-cone effect produces artificial scale-dependence and redshift- 
dependence on the higher-order moments^ of redshift-space clustering of any 
cosmoiogicai objects. 

3. In linear theory we formulated the cosmoiogicai redshift-space distortion^ which 
induces an apparent anisotropy in two-point correlation functions, especially at 
high redshifts. Further detailed studies with N-body simulations E3) indicate 
that it is feasible to constrain the cosmoiogicai parameters from the future 
QSO samples via this effect even though the nonlinear evolution appreciably 
affects the linear theory predictions. 

Apparently, the results described above should be regarded as the first attempts 
to raise the importance and basic features of these two cosmoiogicai effects. These 
results are still far from complete in the sense that there are many aspects which 
remain to be explored. We hope that this short review serves as a practical and 
useful introductory note for more detailed investigations in the future. 
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Fig. 8. Contours of ^ 3 \s±, Sy) for the fixed shape parameter r (= 0.2) and erg = 1.0 models at 
z = 0.2 (upper panels) and z = 2.0 (lower panels); (i? ,Ao) = (1.0,0.0), (1.0,0.9), (0.1,0.0) 
and (0.1,0.9) from left to right. Solid and dashed lines indicate the positive and negative £ , 
respectively. Contour spacings are Z\log 10 |£| = 0.25. 
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Fig. 9. Contours of ^ s \s±, Sy) at z = 3. Upper and lower panels correspond to the COBE- 
normalized CDM models (the shape parameter F is given by Eq. (2-14c)), the fixed r (= 0.2) 
models with as = 1.0. Contour spacings and line types are the same as in Fig. 8. 
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Fig. 10. Evolution of bias and j3{z). Dashed and solid lines correspond to constant biasing param- 



eter b(z) = 2 and 1, respectively, while dotted lines correspond to the bias model (Eq. (2-17)) 
with b — 2 at z = 0. Thick and thin lines correspond to Ao = 0.9 and 0, respectively. 
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Fig. 11. The anisotropy parameter £y ( s )/£± ( s ) 
as a function of z at s — 10/i _1 Mpc in 
cold dark matter universes with Ho = 
70km/sec/Mpc. The upper panels corre- 
spond to a = (empty beams), while the 
lower panels correspond to a = 1 (filled 
beams) . 




Fig. 12. The anisotropy parameter £y 3 ' (s) /£i s ' (s) 
as a function of s at z = 1 (upper 
panels) and at z — 3 (lower panels) in 
cold dark matter universes with Ho ~ 
70km/sec/Mpc. The left panels correspond 
to a = (empty beams), while the right 
panels correspond to a = 1 (filled beams). 
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Fig. 13. £ ' (sx,S||) at z = 2.2 from linear theory (upper panels), and N-body simulations with 
5 x 10 5 randomly sampled particles (middle panels) and with N = 2 x 10 4 most massive halos 
of particles (lower panels). The contour lines represent £' s - ) = 10, 10 ' 5 and 1 (in black), and 
£ (s) = 0.5, 0.3, 0.2 and 0.1 (in white). The white region around the upper-left corner in SCDM 
model indicates < 0. 




Fig. 14. The resulting x 2 -contours in the Oo - Ao plane from the analysis of the data in Fig. 13. 
The theoretical predictions are not corrected for the nonlinear effects. The crosses indicate the 
true values adopted in the simulation models. 
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Fig. 15. P {s) (k s± ,k sll ) at z = 0. The upper and lower panels display the linear theory predic- 
tions and the results from N-body simulations using all particles (N = 256 3 ). The middle 
panels present our nonlinear model predictions on the basis of the non-linear power spectrum 
(Eq. (2-15)). For the values of the pair-wise peculiar velocity dispersion ap, the solid lines cor- 
respond to those from the simulation data, while the dotted lines correspond to the analytical 
fitting formula (3-26). Contour spacings are ZHog 10 P — 0.5, and the contours corresponding to 
10 5 , 10 4 , 10 3 , 10 2 and 10 1 are plotted by the thick lines. 
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Fig. 16. The same as Fig. 15, but at z = 2.2. 
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;. 17. P < - CRD ' , (k a± , fc S ||) at z = 2.2 from N-body simulations using randomly sampled particles of 
5 x 10 5 (lower), 5 x 10 4 (middle), and 5 x 10 3 (upper). The contour levels are the same as those 
in Fig. 15. 
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Fig. 18. The monopole moment, P \ks), of the power spectrum in redshift space at z = 2.2. 
The symbols indicate the mean values of 24 subsamples which randomly select 5 x 10 4 particles 
from the entire set of simulation particles in three different realizations. The error bars are 
computed from the la dispersions among the subsamples for each model. Five theoretical 
predictions are plotted with different types of curves. These values of ap as curves correspond 
to various quoted in the plot. 
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Fig. 19. The constraints on the J?o - op plane derived from the x 2 analysis of simulation data 
at z — 2.2. The upper panels are from P^ CRD ' (k a± , k s \\) with all the simulation particles 
(TV = 256 3 ), while the lower panels are from l CRD \k s ) with 5 x 10 4 randomly sampled 
particles. Thick and thin contours correspond to the results assuming Ao = and Ao = 1 — Oq, 
respectively. The crosses indicate the true values adopted in the simulation models. We assume 
Eq. (3-29) for the value of <t 8 . 




Fig. 20. The same as Fig. 19, but in the fio - us plane. We adopt an analytical fitting formula 
(Eq. (3-26)) for the value of o P . 



